Quantum dark solitons in the ID Bose gas and the superfluid velocity 
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We give explicit connections of quantum one-hole excited states to classical solitons for the one- 
dimensional Bose gas with repulsive short-range interactions. We call the quantum states connected 
to classical solitons the quantum soliton states. We show that the matrix element of the canonical 
field operator between quantum soliton states with A'^ — 1 and A'' particles is given by a dark soliton 
of the Gross-Pitaevskii equation in the weak coupling case. We suggest that the matrix element 
corresponds to the order parameter of BEG in the quantum soliton state. The result should be useful 
in the study of many-body effects in Bose-Einstein condensation and superfluids. For instance, we 
derive the superfluid velocity for a quantum soliton state. 
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The experimental realization of trapped one- 
dimensional atomic gases has provided a new motivation 
in the study of the effects of strong correlations in 
fundamental quantum mechanical systems of interacting 
particles [l|-|3[. Furthermore, localized excitations in 
quantum many-body systems such as in cold atoms and 
optical lattices have recently attracted much interest 
and have been studied extensively in terms of "quantum 
solitons" 0, 0] ■ Localized quantum states are important 
and useful for investigating quite complicated excited 
states of interacting quantum systems. However, it is 
not clear how we can construct or characterize quantum 
states associated with solitons for many-body systems. 
Originally, solitons are special solutions of some classical 
nonlinear partial-differential equations. It is not even 
trivial to see whether there exists a quantum state with 
a soliton-like density profile. 

Let us consider the Gross-Pitaevskii (GP) equation 
which describes Bose-Einstein condensation (BEG) in the 
mean- field approximation Q . We also call it the classical 
nonlinear Schrodinger equation, since it corresponds to 
the classical limit of the quantum nonlinear Schrodinger 
equation satisfied by the canonical Bose field il]{x,t) for 
the one-dimensional Bose gas interacting with the delta- 
function potentials. Here, the system is called the Lieb- 
Liniger (LL) model • The GP equation has dark soliton 
solutions for the repulsive interactions, while it has bright 
soliton solutions for the attractive interactions [1] . It was 
conjectured that dark solitons are identified with Lieb's 
type-II excitations as an excitation branch Q. However, 
it has not been shown how one can construct such quan- 
tum states that are related to solitons or what kind of 
physical quantity can show a property of solitons for some 
states. In fact, each of the type-II eigenstates has a fiat 
density profile, since the Bethe eigenstates are transla- 
tionally invariant. Here we remark that for the attractive 
case, bright solitons are analytic ally derived from some 
quantum states of the LL model [10[ . 



In this Letter we demonstrate that quantum states 
which are tightly connected to classical solitons are con- 
structed from the Bethe eigenvectors of the LL model. 
We call the states the quantum soliton states. Let us 
denote by -(/'qg (x) the matrix element of the field oper- 
ator ipix^t) between two quantum soliton states where 
one state has A^ — 1 particles and another TV particles. 
We show that the matrix element 'i/'qs {x) is well approxi- 
mated by the classical complex scalar field of a dark soli- 
ton of the GP equation in the weak coupling case. We 
suggest that the matrix element ?/'Qg(a;) corresponds to 
the order parameter of BEG in the system with a large 
but finite number of interacting particles in the weak 
coupling case. The result should be fundamental in the 
study of many-body effects in BEG and superfluids. For 
an illustration, we derive the superfluid velocity from the 
phase profile of the matrix element ■0Qs(a;). 

We give remarks. First, superposing Lieb's type II 
excitations, i.e. one-hole excitations, we construct the 
quantum soliton states [ll[, which have broken transla- 
tional symmetry. The Bethe eigenstates are translation- 
ally invariant, while their superpositions are not, in gen- 
eral. Secondly, we show that the amplitude and phase 
profiles of quantum soliton states are consistent with 
those of corresponding solitons of the GP equation. Al- 
though the density profile with a density notch has been 
derived for a quantum soliton state [ll|, the connection 
to solitons has not been shown, yet. Thirdly, it is not 
a priori clear how valid the mean-field approximation is 
for the quantum soliton states. The exact wavefunctions 
given by the Bethe ansatz consist of a large number of 
terms such as A'^!. However, evaluating the matrix el- 
ement ipQa,{x) we identify it as the order parameter of 
BEG. 

Let us consider the Hamiltonian of the LL model [tI: 



N 



1^4 



N 

j<k 



Xk) 



(1) 



2 



Here the periodic boundary conditions (P.B.C.) of the 
system size L are assumed on the wavefunctions. Here- 
after, we consider the repulsive interaction: c > 0. The 
LL model is characterized by a single parameter 7 := c/n, 
where n = N/L is the particle density. We employ 
a system of units with 2m = ft = 1, where m is the 
particle mass. The second-quantized Hamiltonian of the 
LL model is written in terms of the canonical Bose field 
tjj{x,t) as 

^NLs= / dx[^J^^,^^p+c^p^^i>- (2) 

Jo 

where /i is the chemical potential. The Heiscnbcrg equa- 
tion of motion is called the nonlinear Schrodinger equa- 
tion: idtip = —d^ip + 2cip^'ipiil! — fxiil!. 

In the LL model, the Bethe ansatz offers an exact 
eigenstate with an exact energy eigenvalue for a given 
set of quasi- momenta ki, k2, . . . , fcjv satisfying the Bethe 
equations for j = 1,2, . . . , N: 
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Here Ij 's are integers for odd and half-odd integers for 
even A^. We call them the Bethe quantum numbers. The 
total momentum P and energy eigenvalue E arc writ- 
ten in terms of the quasi-momcnta as P 



L l^J = 



E 



Y^^=i ■ If specify a set of Bethe 
quantum numbers Ii < ■■■ < In, the Bethe equations 
([3]) have a unique real solution fci < • • • < fcjy 1121 1 . 
Let us formulate quantum soliton states 
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We shall 

show throughout the Letter that they lead to dark soli- 
tons of the GP equation. In the type II branch, for each 
integer p in the set {0, 1, . . . , iV— 1}, we consider momen- 
tum P = 2'Kp/L and denote by \P,N) the normalized 
Bethe eigenstate of N particles with total momentum 
P. The Bethe quantum numbers of \P,N) are given by 
Ij = -{N + l)/2+j for integers j with I < j < N-p and 
Ij = -{N + l)/2 + j + 1 for j with N -p+1 < j < N. 
For each integer q satisfying 0<q<N— Iwc define 
the coordinate state \X,N) oi X ~ qL / N hy the discrete 
Fourier transformation: 



\X,N) :-- 



1 



cxp{-2nipq/N) \P, N) . (4) 



p=0 



The density profile of the quantum soliton state, 
{X,N\iij'^{x)ip{x)\X,N) versus x, is plotted in Fig. [H It 
is denoted by "Many-body" . Here we have set the coor- 
dinate integer as g = 0, and the density notch is localized 
at X = L/2. The expectation values of the density op- 
erator are effectively calculated [lH by the determinant 
formula for the norms of Bethe cigenstates [l3| and that 
of the form factors of the density operator [14, 
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FIG. 1: (Color online) Density profile of the quantum soliton 
state {X, Nli)"^ {x)rP{x)\X, N) for c = 0.01 and iV = L = 500 is 
shown by a red solid line. The profile of the squared amplitude 
of a dark soliton |'0mf(3;)P with v ~ Vc/2 is plotted with a 
blue broken line. 



The classical complex scalar field of a dark soliton so- 
lution for the GP equation with P.B.C., ^mf{x) (0 < 
X < L), is derived by assuming the traveling- wave solu- 
tion: 'ip{x,t) = iPmf{x — vt). Here we note that the pe- 
riodic soliton solutions of the GP equation are expressed 
in terms of the elliptic integrals We also note that 
the excitation mode has the largest velocity Vc, which we 
call the critical velocity: there is no soliton solution with 
V > Vc 0, H, 13 • In the LL model, the critical momentum 
Pc = mvc corresponds to the Fermi momentum 2k p. 

The density profile of the quantum soliton state and 
the square-amplitude profile of a classical dark soliton 
with P.B.C., ji/'MFla^)!^! agree quite well in the weak cou- 
pling case c <C 1, as shown in Fig. [T] Here we have set 
1; ~ Vc/2 for the dark soliton, and its profile is denoted 
by "Mean-field". 

We suggest that the soliton velocity v ~ Vc/2 is con- 
sistent with the construction ^ of the quantum soli- 
ton state. Each of the type II excitations in the range 
< V < Vc is superposed with equal weight, so that we 
have the average value u ~ Uc/2. It looks like a wave 
packet of the type II excitations. 

Let us consider the matrix element of the field operator 
between the quantum soliton states 

i^Qsix) := {X,N-l\i,{x)\X,N) 

N-lN-2 

(5) 

where P ~ 2ttp/L and P' = 2ttp' /L denote the total 
momenta of the normalized Bethe cigenstates |P, N) and 
\P',N), respectively. We put g = in eq. (O. The 
matrix element {P',N — l\tp{0)\P, N) are evaluated ef- 
fectively by the determinant formula for the norms of 
Bethe cigenstates [l3| and that for the form factors of 



3 



the field operator as 

(P', N - 1|^(0)|F, N) = (_i)A^(w+i)/2+i 



'Af-l N ^ \ I ^ 



n 



det[/(fc,fc') 



v/(fc-,J'+cV v'detG(fc)detG(fc')' 



(6) 



where the quasi-momenta {fci,--- ,fc^r} and 
{fcj, • • • , give the eigenstates |-P, -/V) and 

|P',A^ — 1), respectively. Here we have employed the 
abbreviated symbols kj^e, := kj — ki and k'^ ^ := k'j — k[. 
The matrix G(k) is the Gaudin matrix, whose (j, t)ih 

element is G{k)jj = L + X]m=i ^i^i^m) - K{kj,e) 
for j, ^ = 1, 2, • • • ,N, where the kernel K{k) is defined 
by K{k) = 2c/{k^ + c^). The matrix elements of the 
{N — 1) by (A^ — 1) matrix J7(fc, k') are given by 



Y{aZ,\K-k,+^c) 



nf=i(fca - kj +ic) 



T\a=tj0^a kj) 



{K{k,j) - if(V^)) ■ (7) 



The profiles of the squared amplitude |7/'Qs(a;)p and 
the phase Arg['0Qs(a;)]/7r are plotted in Figs. [2] and |3] for 
TV = 20 and 500, respectively. The squared amplitude 
and phase profiles of periodic dark solitons, iPmf{x), are 
shown by broken blue lines in Figs. [2] and [3] for N ~ 
20 and 500, respectively. They have the velocity v ~ 
2t!-/L. Here the classical complex scalar field tpurix) is 
normalized such that the integral of |'0MF(a;)P with a; 
over the whole region gives the particle number TV. 

The matrix element of the field operator, tpQs[x), and 
the classical dark soliton with P.B.C., V'MF(a;), are in 
good agreement around at the central part of the soli- 
tons in Figs. 2 and 3. In particular, the phase pro- 
files of ipQsix) and ipMFix) completely overlap each other 
(see, the lower panels of Figs. 2 and 3). The profiles of 
square-amplitude |'i/'Qs(2^)P of the quantum soliton states 
are slightly smaller than those of the periodic solitons, 
|V'MF(a;)P, (see, the upper panels of Figs. 2 and 3). For 
= 20 and c = 0.01 the two profiles are proportional 
to each other only with different normalizations. For 
N = 500 and c = 0.01 the two profiles overlap each other 
at the central part and deviate around at the shoulders. 

The agreement of the squared amplitudes |'0Qs(2;)P 
and |'0MF(a;)P should be improved for smaller values of 
c. Let us consider the form-factor expansion of the local 
density at x for the state \X,N): 



{X,N\^p\x)^p{x)\X, N) 

\n)^\X^N-l) 



\{n\ij{x)\X,N)\'. (8) 
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FIG. 2: (Color online) Profiles of the squared amplitude 
i^Qs(a;)P and the phase Arg[i/;Qs(a;)]/7r for the matrix ele- 
ment of the field operator, i/)qs(2;), are shown by red solid 
lines for c = 0.01 and N = L = 20. Those of a dark soliton 
under P.B.C., V'mf(2;), with v ~ 2n/L axe plotted with blue 
dotted lines. Here we put q = Q. 



The second terms of the right hand side of ([5]) give correc- 
tions to \^Qs{x)\^ , and the integral of \'4'qs{x)\^ with x 
over the whole region is smaller than the particle number 
A^. We observed numerically that the correction terms 
become small as the coupling constant c decreases if we 
fix the particle number N , while they increase as TV in- 
creases if c is fixed. In the large TV case, the correction 
terms should be small if the value of c is small enough. 
We thus conclude that the matrix element ■(/'qs {x) is well 
approximated by the periodic dark soliton ipn-pix) in the 
weak coupling case. 

We suggest that the soliton velocity v ~ 27r/L corre- 
sponds to the difference between the average values of 
the total momenta of the quantum soliton states I^T, TV) 
and I X, TV — 1). It is consistent with the structure of the 
matrix element (X, TV - l\'ip{x)\X, TV). 

We now argue for the claim that the matrix element 
V'QS {x) gives the order parameter of BEG in the quan- 
tum soliton state \X, TV) for the weak coupling and large- 
TV case. Here the system size L is also very large since 
wc set n = N/L = 1. We denote by pi(a;,?/)|^) the 
onc-particlc reduced density matrix for a given state 
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FIG. 3: (Color online) Same plots as in Fig. [2]with the system 
size N ^ L = 500. 



I^*): pi(a;, y)|5,^ = [x)Tip{y)\^) . Wc now conjec- 

ture that the matrix element ^qs{x) satisfies the rela- 
tion: pi{x,y)\x,N) - '0Qs(2;)^Qs(2/) for |a; - y| > 1. We 
call it conjecture A. Here we assume that the system size 
L is much larger than the healing length Ic = 1/ \fcn. For 
instance, Ic = 10 for c = 0.01 and n = 1. 

Suppose that the one-particle reduced density matrix 
for a given state j^") is diagonalized as 



pi{x, y)\^) 



(9) 



For the ground state we can numerically show that for 
small c the largest eigenvalue uq of pi {x, y) is much larger 
than the other eigenvalues: uq ^ rij for i 7^ 0, i.e. the 
existence of BEC ll|. For the state \X,N) it could be 



technically nontrivial to diagonalize pi{x, 'y)\x,N) numer- 
ically. However, instead of doing it we point out that 
conjecture A is consistent with the following observation: 
The quantum soliton state \X, N — 1) is dominant among 
the intermediate states in the expansion ([5]). Here, we 
estimate the fraction of the correction term from the dif- 
ference between the local density at x for the state \X, N) 
and the squared amplitude |-0qs(x)P, and it is small for 
small c and large N. Moreover, from the difference we 
estimate the condensate depletion, i.e. the fraction of the 
non-condensate components. It should have the largest 



values ioT X = y, since the local density at x gives the 
diagonal element of pi{x, y)\x,N) with x — y. 

We therefore conjecture that for small c and large 
the order parameter ^/n^Xoi^) of the quantum soli- 
ton state \X,N) is given by the matrix element ipQs{x), 
which is well approximated by the periodic dark soli- 
ton ipMpix) with V ~ 2'k/L. Here, the order parameter 
y/noxo{x) has been defined by eq. (|9]) for j^f) = |A, A^). 
In terms of BEC we have connected the dark soliton 
iPmf{x) with V ~ 2tt/L to the state |A, A^). It was not 
trivial to specify the soliton velocity v. 

Let us now derive the superfluid velocity for a quantum 
soliton state |A, A^). For large A^ such as A^ = 500, the 
phase field is fitted by d{x) = ttx/L — ttH{x — X — L/2), 
as shown in the lower panel of Fig. 3. Here H{x) de- 
notes Heaviside's step function: H{x) = 1 for x > 0, and 
H{x) = otherwise. Numerically we observed that the 
phase profile does not depend on the value of c for large 
A and small c. We derive the superfiuid velocity from 
the phase field of the macroscopic wavefunction, 9{x), by 
Vs = 2{h/2m){de/dx) For large A^ we thus have 



2tt 

— — -K 5{x 
Ij 



X- L/2). 



(10) 



The superfluid velocity Vs has a singularity at the lo- 
cation of the soliton, and is consistent with the soliton 
velocity v ~ 2tt/L. 

The finding in the Letter suggests several possible fu- 
ture researches in quantum dynamics such as the collision 
of two quantum solitons, which is nontrivial in the dy- 
namics of interacting quantum systems [lij . 

In conclusion, in order to prove that the quantum 
states |A, A^) constructed in (j4]) for the ID Bose gas are 
closely connected to classical solitons, we have shown the 
following two points: First, the density profile of the state 
I A, A^) is consistent with the profile of the squared am- 
plitude |'0MF(a;)P of the periodic dark soliton of the GP 
equation with v ~ Vc/2. Then, the matrix element of 
the Bose field operator, 'ipQs{x) (A, A^ - l|-0(a:)|A, A^), 
coincides with the classical complex scalar field V'MF(a;) 
of the dark soliton of the GP equation under P.B.C. with 
V ~ 2tt/L. The agreement is good for small c. Further- 
more, we suggest that the matrix element ipcis{x) gives 
the orer parameter of BEG in the quantum soliton state 
|A, A^) for small c and large A^. 

The authors thank I. Danshita and K. Sakai for their 
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